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1.0 SUMMARY 


II I' ultimate objective of work identified an Tusk I of NAS2*«729 is to develop and 
implement, within the FLEXSTAB computer program system, improved aerodynamic 
mode I a and methods to replace those currently embodied in the FLEXSTAB aerodynamic 
influence coefficient 'Alt’' program 

The work is to be based upon the aerodwiumic technology methods and experience which 
have evolved from the vortex spline scheme developed originally under NASA- Ames 
contract NAS2 d.130 and the subsequent contract work w hich has gone beyond the original 
vortex spline to the doublet and source spline concepts, work presently being conducted 
unde. NASA-Atnes contract NA82-7729 ( Development of a FLEXSTAB Computer 
I’ropran It is the latter work, in particular, which supplies the impetus, encouragement 
and strong technological nasis to enter into the Task I work leading toward improved 
uerodynamic methods for the FLEXSTAB computer program system. 

The conclusion drawn from this studv is that the ultimate objectives of Task I can he 
achieved with little technical risk and without loss of any capability relative to the 
NASA-Ames released version I 01.00. The direction of that development and the ability 
of the FLEXSTAB computer program -vstem to interact dynamically with advancing 
technology hinges upon development of a restructured program system. Two options are 
presented 

Option A assumes that the present FLF.XSTAB system will not undergo a program 
restructuring 

Option B assumes that FLEXSTAB will be restructured. The exact body surface 
aerodynamic model is recommended to provide an improved modeling capability not 
presently available. 

Option B is the specific recommendation of this report. 


2.0 INTRODUCTION 


Steady state aerodynamic technology has advanced beyond the original numerical 
littmK surface technology of con slant pressure parti- I k introduced by Woodward (refs I 
through 3>; ins aerodynamic model for wing-body interference ( interference shi-ll of 
constant cross section) has also been superseded It is the Woodward model which is 
employed in the current KLKXSTAH program system. New technology, however, has 
been and is being developed which has the potential to improve the aerodynamic 
methods within KLRXSTAB Chief among the characteristics of the new technology 
are the capability to treat theoretical models that more closely represent the true 
configuration surfaces, insensitivity to configuration paneling (and the associated 
promise of fully automatic paneling), computational simplicity (computation of Alt' 
integral expressions in closed form), and the possibility of producing potential flow 
information that will provide an adequate foundation for future implementation of drug 
force prediction involving houndary layer analysis (drag prediction in KLKXSTAH is 
limited to induced drag'. 

The work described in this report is an analysis of alternate approaches which presents 
the relative merits of each of several possible candidate theoretical aerodynamic models 
and numerical schemes, now or soon to be available Kach of the schemes is evaluated 
in each of the following specific areas: (1) applicability to the steady flow needs of 
KLKXSTAH (subsonic and supersonic) and (2) adaptability to (u I the low reduced 
frequency approximation, * h i Kussner and Wagner functions approximation, and (c) 
arbitrary reduced frequency analysis. Existing and new technology is defined and 
reviewed. 

The relative merits of each of the various methods are assessed with respect to their 
features and limitations as well as their practical applicability to the present and the 
potential growth versions of KLKXSTAH. Particular attention is directed toward the 
steady flow and low reduced frequency requirements of KLKXSTAH, with a view toward 
adaptability to arbitrary reduced frequency analysis. Consideration is given to the 
potential gains in accuracy measured against the computational cost of obt r a 
solution. 



3.0 SYMBOLS AND AHiil. KVIATIONS 


A1C aerodynamic influence coei jnt 

a t Npeed of sound 

centerline 

C,, dimensionless pressure coefficient 

I) IH total derivative 

dx'.dy' differential elements 

ES1C external structural influence coefficient 

F equation of bounding surface 

KiSi singularity strength 

ISIC internal structural influence coefficient 

K kernel function form 

k ratio of specific heuts 

LI) line doublet 

LS line source 

M Mach number 

f\P x static pressure at arbitrary point and at infinity 

magnitute of total velocity 
t time 

U x magnitude of fluid velocity at infinity 

VC constant vortex 

x,y,z coordinate axes ix. streamwise direction) 

V gradient operator 

d> perturbation velocity potential function 

<t> total velocity potential function 

/*x fluid density at infinity 

iT,n,y singularity strength distribution function for source, doublet, and vortex 

r) time derivative 

in 


t 
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4.0 OVKKVIKW OF FLKXSTAH 


4.1 SOME PRESENT FEATURES* 


(•oner a I Description 

FLEX8TAB is a collection of computer program* ••hiiii have been assembled for 
aeroelastic analysis of arbitrary airplune-like configurations The programs can he 
linked by tape or disk data transfer. These programs can be executed singly or linked 
for consecutive execution in a single run Thus, a complete ueroelastic analysis made up 
of an aerodynamic analysis, structural analysis, calculation of static and dynamic 
stability derivatives, time history calculation, etc., can be accomplished with minimum 
user participation leaving the complicated interface calculations to FLEXSTAB 
(see fig 1). 



Boundary Value Problem 

The aeroelastic formulation of FLEXSTAB is based upon linearized partial differential 
flow equations, the associated How boundary conditions, and a linearized structural 
model of the configuration (nonlinear effects can be introduced at the problem analysis 
level, e g., the introduction of wind tunnel derived-s'ahility derivatives into the time 
historv calculation). A complete description of the aerodynamic boundary value problem 

rinoco. K N . and Mercer, .1 K. FLEXSTAB - A Summary of the Functions and Capabilities of 
the S'ASA Flexible Air/dane Analysis Computer System Hoeing document D6-41098, NASA 
CR-2564, December 1973 
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which KI.KXSTAH addresses is contained im section fi.O. The aerodynamic and 
structural models are based upon deformation* producing linear presMures, and forces 
producing linear deformations, respectively. Great simplification results and arbitrary 
motions can be composed from superpositions of linear solutions 

Numerical Scheme Aerodynamic Modeling 

The numerical models currently employed in '"'LKXSTAB are nuinulaniv solutions of 
the partial differential How equation The> .ire >f three has’'' forms within KI.KXSTAH 
source, doublet vortex The source solution is used to represent thickness eflects; the 
doublet and vortex to represent Idling effects In general, the thickness Htid lift effects 
can be separated because of tin* assumed linearity of the problem These singularity 
solutions are distributed over regions of space and are constructed as line segments or 
panels The solutions automatically satisfy the differential equations. A complete 
analysis is accomplished bv linearly superimposing combinations of the singularity 
distributions with the strength of each distribution being determined from the boundary 
conditions at specific boundary points 

Body like elements are composed of u line singularity distribution >f sources to 
represent the thickness, and a line singularity distribution of doublet 4 to represent 
lifting effects The source distributions are constants over segments of the line in 
subsonir flow and linear in supersonic flow, the doublets an- quadratic functions 

Wing-like elements composed of panel singularity distributions of sources 
(thickness) and vortex elements (lifting effects 1 Tin- sources are linearly varying in the 
direction of the free stream; vortex distributions are constant over each panel 

Wing-body interference (lifting) effects on the bodies are carried by cylindrical shells ol 
constant cross sections which shield the hody-like elements from all other elements. 
Vortex panels are distributed over these shells, each pr.nel has a constant strength On 
the wing, constant vortex panels carry both the interference and lifting effects 

The leading-edge suction at thin wing leading edges must be computed to provide their 
contribution to the induced drag (negative contribution). This is done in KI.KXSTAH by 
application of the principle of conserved momentum. 

Wing-body intersection is treated in a straight forward manner. The wing i vrcepts the 
body ut the interference shell and the intersection line always occurs ai ? inel edges. All 
carryover interference effects are provided by the vortex panels on the shell and wing 
No special treatment is required to provide thi appropriate jump in pressure across the 
wing surface at the intersection because <T the constant vortex (or pressure) panel 
elements 

Nacelles are treated in several ways In each case they are modeled with line 
distributions of sources and or doublets and with an interference shell which may (user 
option) experience interference effects. 

Struts are treated as wing-like surfaces Struts connect nacelles to either wings or 
bodies. 
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Limitations «»r Restrictions of the Numerical Scheme Aerodynamic Modeling 

I'ln partial dilTerential flow equations und How boundary conditions, along with the 
Ntruetural model arc linearised This can be a serious restriction for some flow 
conditions. The linearised wing assumptions break down for thick wings and 
interference effects would accordingly be incorrect. Panels ar required to have their 
side edges parallel to Ihe direction of the I’ > • stream and the vortex panels have an 
inherent planar wake extending to downs. ream infinity, the edges of which have 
discontinuous I infinite! jumps in perturbation downwush velocity. This imposes 
restrictions < n panel size and spacing, with downstream surfaces dictating paneling 
requirements on upstream s irfaces. Boundary point placement can be crucial A dense 
control surface paneling may require an overly large nuintier ot panels over the 
remainder of the configuration 

Ihe interference shell represents a mean surface ol the body Accordingly boundary 
conditions are not satisfied on the actual body geometry Results may be adequate for 
bodies of nearly circular cross sections hut can introduce serious error for more general 

bodies. 

It is often risky to place small vortex panels in proximity to large ones. Near a 
wing-strut intersection this may be unavoidable and nunn rical errors can be introduced 
in these regions. This can also occur when paneling the wing in the spanwise direction 
(stream direction panel spacing is not as critical). 

Computation 

Very large cases hav been run in FLEX8TAB. which is virtually an open ended 
program. Economic considerations, however, tend to keep the number of aerodynamic 
singularities to n reasonable level (• 500). It is not entirely realistic to quote 
computation times for the aerodynamic influence coefficient (Alt') as they depend on 
many factors Different schemes can only be compared when referenced to exactly the 
same configuration. Such times are therefore not quoted here 


4.2 FUTURE GROWTH POTENTIAL 


Ultimate Objectives 

The ultimate objective of Task I is to improve portions of the subsonic und supersonic 
aerodynamic programs by replacing them with a new FLEXSTAB Alt' program 
consistent with the theoretical aemelastic formulation currently in FLEXSTAB (but 
founded upon an advanced generation ol aerodynamic building blocks) and to implement 
these advanced aerodynamic methods into FLEXSTAB. The advanced methods will 
evolve from adaptations of technological developments conducted under NASA-Ames 
contracts N AS2-H5.IO and NAS2-7729. The end product will be a new version of 
FLEXSTAB having improved aerodynamic features that provide greater accuracy and 
computational efficiency; and that allow more flexibility in the size, shape, and 
arrangement of aerodynamic panels. This new version of FLEXSTAB is to have no loss 


b 



m 'inv capability relative to the initial NASA Amen rrlritwi version (1.0100) In 
particular a low frequency type* approach to the unsteady *«*•■< dynamics is to be 
retained 

l.imitwtinns to FI KXS'I All Growth 

A* presently constructed. FI.KXSTAH ih not a system of independent computational 
modulc'N; it ih rather a closely coupled system wherein a change in one portion affect h 
the others. For example, the constum vortex panel Alt’* could not be immediately 
replaced with another numerical model without serious consideration of the impact on 
the remainder of FI.KXSTAH For this reason. two optioiiH are being considered for a 
new version ol FI.KXS TAH The lust. designated option A ."lams the present structure 
and aerndvnumic modeling of FI.KXSTAH (line body, mean surface wing, and 
interference shell!; option H involves the restructuring of FI.KXSTAH and use ol the 
exact body representation. The two options will be discussed in section 8.0 
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5.0 BOUNDARY VAI UK I’KOBI.KM INSCRIPTION 


VI »*Aii . IAI. llimUKM I Al. EQUATION 


The par. ml differential flow equation (ref 4* in derived Crum the Kul«*rmn momentum 
f , lation, continuity, and a r»'li»tinn fur the speed ol Mound undor the HNNumpiionN of 
perfect fluid, no external force field, no heat conductivity, continuum, inviacid, 
adiabatic, irrotational, except for hounding surfaces and certain prescribed regions 
T‘>ih equation In linearized assuming rertain orders of magnitude for the perturbation 
velocity component** The resulting partial differential flow quatmn, with axes fixed to 
the v infiguration. in 


where 


*^\X * 1 M * I 4 lUyy 4 <*> 


2N1 I 

X* ' *Xl "j <^tt 

• x m% 


III 


«x 

t 

M 


coordinate uxch ix, atreamwiaa direction) 
|>erturbation velocity potential function 
speed of Nound 
time 

Much number 


Equation (I) in the governing equation for the aerodynamic influence coefficient 
program within FLFXSTAB as well no for the advanced aerodynamic building blocks 
Thi* atatic preaaure at a point in the fluid can be related to the free atream Htatic 
pressure anil dynamic pressure by defining the dimensionless pressure coefficient, Up, 


where 


c. 


p p» 

i' 


P. P* 

f*x 

Ux 


Htatic pressure at arbitrary point and at infinity 

fluid density at infinity 

magnitude of fluid velocity ut infinity 




An exact expression for C„, derived under the assumptions listed, written in terms of 
the velocity potential function is I ref. 5): 


k 



where 

k ratio of specific heats 

<l> total velocity potential function 

Q magnitude of total velocity 
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I hi* I'Uil Op expression I<»r an incompressible fluid i*» 



Each *•!' the expressions tor the pressure coefficient t(3>, <40 are independent of the 
ci-ord nate system used These equations can he linearised, in which case the resultant 
expressions may not be independent of the coordinate system. In addition the linearised 
expressions will differ depending upon the assumed orders of magnitude of the 
perturbation velocity components FLKXSTAB uses a linear expression for <‘ p 

5.2 HOLNDAKY CONDITIONS 

It F <x, y, /. ti 0, i the equation of the bounding (configuration i surface, then the 
general boundary condition is < ref ti>: 

. o - — ♦ V* • VF <5. 

1)1 dt 

and is tin relation winch limits the flow to he tangential to the hounuing surface The 
boundary condition can he linearised and various expressions obtained depending upon 
the assumed orders of magnitude of the perturbation velocity components and the 
ciMi-dmate system Kl.KXSTAB uses linearised boundary conditions 

5.3 NUMERICAL SCHEMES 

Many schemes and techniques have been proposed to solve the aerodynamic boundary 
value problem The most promising of the methods for aoroelustic analysis are those 
which supply sufficient information for the structural representation without 
unnecessary comple xity or over simplification. The methods having the most promise for 
the Kl.KXSTAB uemelastic analysis are smKulanty methods For completeness, 
additional methods which have been used for the aerodynamic boundary value problem 
are listed below us other methods 

Singularity Methods 

Singularity methods have their foundation in the existence of fundamental solutions of 
the linear flow equation <1) (ref. 71. Derivatives of fundamental solutions and linear 
combinations of them ulso satisfy equation (1). From this knowledge are derived the 
uerodynumic singularities: source, doublet, and vortex These singularities un- 
distributed on regions oi the bounding surface; their value or strength is determined 
from tin- boundary conditions on the bounding surface. 

The singularity methods have taken various forms. They have been applied to steady 
and unsteady and subsonic and supersonic flows <refs H and 91. They have been placed 
at points along lines, und over area regions of the bounding surface I refs. 1, 10, and lit 
Their strength variations over regions of the bounding surface have been constant, 
linenr. and quadratic (refs 12 through 14>. They have been made to span single and 


muUir.li- line segments and single and multiple urea regions ias well ss the entire 
eonf.guration < rein 14 through Hi). The recommended forms of the singularity solution* 
under consideration lor improving the FLEX STAB Alt’s will be presented in 
.•action 7.0. 

Other Methods 

Among the other methods which have been used for aerodynamic boundary value 
problems are ire*s. fi and 17 through 22 ): 

Finite difference 
Finite element 
2-, :i-l) characteristics 
Strip theory 
Slender body theory 
Conical flow theory 
Newtonian impact theory 
Piston theory 

None of these other methods are considered adequate for a complete FLEXSTAB 
aeroelastic analysis 


5.4 AERODYNAMIC MODELING 

The following considerations are necessary when modeling an aerodynamic configuration 

Body 

It is highly desirable to model the actual body surface, something which FLEXSTAB 
does not do. The interference shell would no longer be necessary and u more realistic 
discretisation of the bisiy would be obtained in terms of surface boundury conditions and 
resultant solution | fissure distribution!. This modeling has been achieved with success 
i ref. 2o> and has been applied to a wide variety of configurations. 

Wing 

The wing can be modeled ns a thick surface (actual surface discretization) or chin 
surface (surfaces across which the surface velocities are approximately parallel). Ii. 
subsonic flow, each of these models has been applied successfully. In supersonic flow the 
hin wir.g approximation is considered appropriate, because of the existence of 
physically thin surl-vrcs and numerical difficulties anticipated when dealing with the 
exact surface.-, in ton-close proximity in supersonic flow. 

Leading- E« ige Corrections 

Lcadinx-edgc suction is the term designating the component of induced drag which 
cannot be obtained by resolution of forces over a thin wing surface in potential flow. It 
can be obtained bv application of the principle of conserved momentum It arises from 
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the requirement that the How turn through a finite angle in an infinitesimal distance 
The result is an infinite pressure at the (subsonic) edge In turning, the infinite 
pressure acting over an infinitesimal area produce* a finite contribution to the induced 
drag, u component which actually exists on a real configuration for which the leading 
edge ha* a finite radius 

A thin aing having a sharp leading edge may not have fully attached How on the edge 
The leading edge auction analogy of Polhamus (ref. ‘24 1 can be used to obtain the 
incremental nonlinear lift resulting from the separated vortex flow. 

A subsonic leading-edge correction has been formulated iref 25) mr a thick wing with 
linearized boundary conditions. 

Wmg-ltndy Combination 

It is often difficult to model the geometry of a region where a wing and body come 
together Thin wing assumptions may break down in these regions and local surface 
pressures may be in error. A program in which the body can be accurately modeled 
helps eliminate these problems In addition, the treatment of the carryover lift must he 
properly modeled 

Nacelle 

It is important to correctly model the inlet and exit Hows of a nacelle. In order to 
accurately predict pressures near a nacelle and the interference effects which it creates, 
an adequate geometric representation and numerical scheme must be available. 

Strut 

Struts are usually small in size and may be located relatively close to regions of large 
numerical influence It may be difficult to model them accurately. 

Wake/Exhaust 

The accurate modeling of wakes and jet exhausts may lx* required in order to represent 
an important feature of the physical flow as, for example, the wake from a high lift 
surface The proximity of the wake or exhaust to control surfaces is important for the 
calculation of control deiivativea, tail sizing, etc. Although the wake surface position is 
generally unknown, a more accurate placement of the modeled wake might be available 
by iterat ion of this nonlinear problem, experimental data, or intuition. 

Control Surface 

The numerical treatment of control surfaces is important to obtain accurate control 
surface characteristics. Depending upon the numerical method selected, the control 
surface modeling may be straightforward or require special model tailoring to avoid 
numerical instabilities. 


0.0 CANDIDATE AERODYNAMIC MODELING SCHEMES 


Several aspects of aerodynamic modeling were presented in section 6 4 F.ach of the 
items has important bearing on the full boundary value problem soli.lon and, therefore, 
is given special attention. A discussion is presented here jf the specific areas of 
aerodynamic modeling, the types of models, and their respective features. 

Body 

Two of the methods to model the body are ill the linearized boundary condition 
interference shell concept, coupled with separating the body into the slender body lifting 
and nonlifting problems t introduced bv Woodward (refs. 1 through 3) and subsequently 
used in KLEXSTAB (ret 26>). and (2) the exact boundary condition surface paneling 
technique ref 23). These methods have not only been used to model fuselage-like bodies 
but also nacelles. Nacelles are discussed under a separate heading 

Consider first the linearized problem The actual body may not be slender nor have a 
circular cross section. However, these approximations may still be made. The 
approximation is often made that the (noncircular) body has a circular cross section of 
equi valent area In addition, the slender body approximation allows splitting the flow 
problem into the two linearized problems of axial How over a body of revolution (the 
nonlifting problem) and u cross (low over the same body u..e lifting problem). For 
slender isolated bodies of revolution in potential How, the results are good. Woodward 
(refs 1 through 3). and Dusto (ref. 2t>( have both coupled this approximation with a 
cylindrical interference shell to perform wing-body analysis. The linearized 
approximation has application to both subsonic and supersonic How The results can be 
quite good for specialized wing body combinations but a definitive statement regarding 
a general configuration cannot be made 

Alternatively, in subsonic flow, the actual body surface can be discretized into surface 
p nels to provide a more accurate description of the physical geometry. The governing 
flow equation and the surface boundary conditions are linear. Paneling the body surface 
in supersonic Hows is still not generally available, although Woodward reports 
successful results (ref. 27). 

The coupling of the linear aeroelastic analysis and the linear aerodynamic problem is 
through the aerodynamic load. The use of a paneled body surface in the aeroelastic 
analysis will require additional study because the aerodynamic boundary value problem 
is linear in the potential function and the perturbation velocities, while the load is a 
nonlinear function of the velocities if the load expression has not been linearized. 

Note: For reference purposes, the following will be usef in section 8.0 (see fig 2). 


i: 


U.1 

Line representation 



8 2 

Surface representation 



Figure 2. Notation for Aerodynamic Models of Body 


Wing 

The wing may be modeled in one of two ways. First, by linearizing the boundary 
conditions, the wing geometry and boundary conditions muy be represented in some 
mean surface as u decomposition into the nonlifting part (thickness! und lifting purl 
(camber). This is a formulation which lends itself directly to the linear aeroelastic 
problem The reaulting wing-alone solution for a thin wing shows agreement w ith more 
exact methods both in subsonic and supersonic flow, with the possible exception of the 
subsonic leading edge For a thick wing, the linearized solution may still be useful 
although surface pressure may be expected to be in error in some regions. This is 
general I. \ unimportant for aeroelastic analysis where the aerodynamic load is of 
primary importance. 

Secondly, the actual wing surface may he paneled and exact boundary conditions 
applied at the wing surface For thin surfaces, there are two objections to this approach 
(1) the solution may be numerically unstable for closely spaced surfaces whose surface 
perturbation velocities are approximately parallel and <2> the number of panels is 
doubled, rr is the size of the wing influence coefficient matrix For thick wings, the 
surface paneling approach may he the only way to achieve satisfactory accuracy for 
detailed surface pressures Ir. subsonic flow, this technique has proven successful It has 
not yet been successful in supersonic How. although it has been investigated to a 
limited extent (ref 27). 

Nolo: For reference purposes, the following will be useful in section 8.0 (see fig 3). 


W.1 

Mean surface representation 


W.2 

Exact surface representation 



Figure 3. Notation for Aerodynamic Models of Wing 


Leading Kdge Corrections 

If the wing is analyzed as a thick surface, by using surface paneling and exact surface 
boundary conditions, there is no need to separately calculate the leading edge suction 
(or thrusti. However, as mentioned in section 5 .s' , for a thin surface and linearized 
boundary conditions, it must be calculated separately to obtrin an accurate value for 
the induced drag. For this case, it may also be necessary to introduce a calculation to 
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obtain accurate value* of surface pressure near the subsonic leading edge to properly 
account lor the thickness effect* (e.g., kit-gels' rule, (ref 25i). 

The simplest method to obtain the leading-edge suction ih to use the conservation of 
momentum principle. It requires perturbation velocities in the vicinity of the subsonic 
leading edge, a string, urward calculation with any of the numerical schemes 
presented in section 7.0 (for verification of application see ref. 24). 

Wing-Body Combination 

The isolated body and wing model formulations were discussed in the foregoing under 
their respective headings Any combination of those body and wing aerodynamic 
modeling schemes could be combined for th<- wing-body problem in u FLEXSTAB type of 
environment Tin combinations which have been successfully demonstrated for subsonic 
and supersonic flow are showm in table l (refs 1 through 3, 26, and 27). 


Table 1 .— Wing-Body Aerodynamic Modeling Schemes 


Body 

Wing 

Leading edge 
corrections 

Comment 

Line * 

shell 

representation 

Exact 

surface 

Mean 

surface 

Exact 

surface 

X 


X 


X 

Subsonic 

Supersonic 


X 

X 


X 

Subsonic 

Supersonic 


X 


X 


Subsonic 


Lawrence and Flax ref. 28i report some of the various formulations which have been 
used to represent wing-body carry-over effects. These methods are often semi-empirical 
due to the difficult nature of the problem. (The wing is assumed to be attached to a 
cylindrical body of constant cross section and infinite length.) The methods are 
specifically tailored to the flow region, the aspect ratio, and other geometric constraint 
Except for the integral representations of the wing-body problem, none of the other 
formulations presented in the foregoing reference are considered to be candidate 
aerodynamic modeling schemes in a FLEXSTAB type of environment. 

For FLEXSTAB application, the following modeling schemes for the carry-over are 
being considered. The first is the current interference shell representation of Woodward. 
In this model, a mean cylindrical shell of constant cross section is used to model the 
interference effects. The body ia line singularity representation! is considered to be an 
isolated body in a uniform freestream. 

The second scheme is an interference shell of nonconstant cross section. 

The third scheme, not properly termed an interference system, is an internal lifting 
system that is an extension of the wing into the body. 
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Note Kor reference purposes, the follow ui)< will be useful in section 8.0 (see fig 4) 


1.1 

Interference shell of constant cross section 
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1.2 

Interference shell of nonconstant cross section 
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Internal liftinq suiface 



Figure 4 Notation for Aerodynamic Models for Interaction Effects 


Nacelle 

The (low characteristics in the region of a nacelle (ref. 29) are probably the most 
difficult to predict and will be strongly influenced by the aerodynamic modeling scheme 
selected. A nacelle is generally ir close proximity to other surfaces (wing, body, strut, or 
another nacelle) ao that the aerodynamic modeling of the various surfaces cannot be 
done independently 

There are certain "haracteristics of nacelles which should be considered, among which 
are the following: 

• Inlet characteristics are important and often must be preserved to guarantee 
powerplant performance. These include: proper inlet velocity ratios and inlet 
streamline directions. 

• At the exit nozzle it may be necessary to apply the Kutta condition. 

• Leading-edge (inlet lip) corrections may be important. 

• Exit velocity ratios and streamline directions affect interference pressures. 

• Exhaust entrainment and displacement characteristics are particularly important 
at low speed. 

• Momentum changes of the fluid passing through the engine produce forces on the 
nacelle w hich cannot he- calculated by resolution of surface pressures. 

• It may he necessary to properly account for external geometry discontinuities 
associated with exotic nacelle designs 

• Inlet unstart simulation is very important to the stability and control 
characteristics (yaw and roll) at high supersonic Mach number ( -2). 
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For FI.KXSTAB application, the following modeling schemes arc considered candidates 
for nacelle molding The firm ih a simple How through nacelle The nacelle in modeled 
as an open body; that ih. a body with a hole in it It may have either a constant or 
nonconstant cross section. No means ih provided to control the inlet or exit velocities, 
they are dictated by the Kutta condition applied at the exit 

A nacelle can also be modeled hy representing it uk a stream tube whose geometry ih 
specified and which ih allowed to deform in an upp'oprate manner This ih termed a 
streamtube nacelle 

The third model is a simulation of an engine-on condition. This model has a barrier 
across the (open) nacelle on which the muss (low is specified. 

A schematic of each of these is shown in figure 5, along with notation used in 
section 8.0. 
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Figure 5. Notation fur Aerodynamic Models of Nacelle 


Strut 

A strut is generally small in size and located close to other surfaces (wing body, or 
nacelle) It may be thin or thick and have either linearized or exact boundary 
conditions, as was discussed under the wing heading Its physical shape may strongly 
influence the flow characteristics and its design may be used to control the nearby flow 
characteristics. 

Wake/Exhauet 

One wake model and one exhaust model are considered candidates for FLEXSTAB 
application The wake model, originating at lifting surface trailing lor leading) edges is 
specified; that is. its position is assumed to be known The wake is defined by the 
geometry of the trailing surface. The plunar wake is a subset of thiB more general 
model. 

The exhaust model is intended to simulate the existence of a plume emanating from a 
nacelle. A plume shape and entrainment distribution are assumed to be known 

The wake and exhaust models are illustrated in figure 6. along with the notation for 
these models which is used in section 8.0. 
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Figure 6. Notation for Aerodynamic Models of Wake and Jet Exhaust 


Control Surface 

The control surface# can be very adequately modeled with the previously mentioned 
modeling techniques lor the wing and wake. In addition to all of the previous 
discussions, it is important to have the option to model a complete configuration, as 
contrasted with modeling schemes leg, Fl.KXSTABi which assume certain «vmmetr> 
properties for the configuration and deul only with a half model. Such an option would 
enable analysis of antisymmetric configurations such as the slewed wing iref. 30). 
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7.0 CANDIDATE NUMERICAL SCHEMES 


7.1 SELECTION Ol CANDIDATE SC HEMES 

The numerical scheme* considered candidates for FLEXSTAR all full into the class 
designated singularity methods The fundamental solutions of the linearised flow 
equation ( 1 ), are 


Source- 

</» (x.y.z.t) 

Si- 

s' 

(x’.y'.ti K <f (x',y'; x.y.z.ti dx'dy' 

(6) 

Doublet: 

ib (x.y.z.t) 

II M 

S' 

(x .y'.ti K^'x’.y'; x.y.z.t > dx'dy' 

(7) 

Vortex 

•b (x.y.z.t) 

Ur 

S' 

ix'.y'.ti Kyix'.y’; x.y.z.t* dx'dy 

(Hi 


Each solution <b (x. y. z, and t) is expressed us an ii tegrul relation involving an 
unknown singularity strength distribution function (<r. p, and yi und a known kernel 
function form iK, f , K^. und Kyi The singularity distributions are shown written for the 
plane Z' 0 The kernel functions depend upon the flow regime (subsonic or supersonic) 
and whether the flow is steady or nonsteady. However, they are not shown here since 
the exact forms do not contribute* to the discussion, but may be found in references 4 
and 31 . Directly of interest to this discussion are the forms of the strength distributions: 
it, m. and v. 

Source 

Source singularity solutions have probably been used for more varied applications than 
uny other singularity method. This is due first to its mathematical simplicity and 
second to its ability to simulate a real, physical flow They have been used to model 
complete wing-body combinations (ref. 23), nonlifting subsonic flow (ref 32) and lifting 
supersonic flow (ref. 33 1 

A surface distribution of sources has a jump in normal velocity associated with crossinn 
from one side of the surface to the other. Sources have been used to represent notwuling 
thick bodies or wings for both steady and nonsteady flows. Supersonic lifting wing 
problems can also be formulated with sources such that the upper and lower surfaces do 
not experience one another's influence (Mach box (refs. 33 through 36) and 
characteristic box (ref. 36)). 

In its simplest form, the source strength (<r) is assumed to he constant on a given region 
of the bounding surface. In general, the strength can he taken arbitrarily, having the 
parameters which define the variation determined in a specified manner by conditions 
on the hounding surface. 

As a specific example, the source spline scheme is being formulated and developed 
under Task III of NASA-Ames contract NAS2-7729. It has two characteristics. First, the 
strength variation is linear in each of two coordinate directions. Second, the surface 
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curvature effects arc included within the rurniulatnm The constant strength, Hat panel 
representation ( refs 23 and 37) ih available us a special case The practical application 
of the source spline is easily recognized and the consistency of the mathematical 
formulation hus been demonstrated (ref. 38). It has application both to subsonic and 
supersonic flow, both to lineatized and exact boundary conditions 

Doublet 

A surface distribution of doublets has a jump in tangential velocity associated with 
crossing from one side of the surface to the other Such a singularity surface finds 
application in lifting problems where tumps in velocity potential occur, e.g.; subsonic 
and supersonic, and for both steady and nonsteady The doublet is a natural numerical 
modeling element for the wake region associated with lifting solutions 

The doublet strength ijii may be taken constant (Doublet Lattice, refs 39 and 40) or 
have some nonconstant variation (doublet splines Tusk III. NAS2-7729). Although 
mathematically more simple, the doublet lattice is only a subset of more general 
distributions, such as the doublet spline 

As a specific example, the doublet spline is being formulated and developed under Task 
III of NASA-Ames contract NAS2-7729. It has u strength variation which is quadratic 
in each of two coordinate directions and includes surface curvature effects within the 
formulation The vortex lattiei type of numerical method is available as a special case 
No other reference to the doublet spine is known. 

Vortex 

Vortex models have an inherent planar (kinematic) wake They are formulated by 
superimposing doublet solutions and are useful for lifting problems where the physical 
location of the wake is relatively unimportant, as is the case for the linearized, thin 
wing problem The surface vorti-x distribution has a jump in tangential velocity 
associated with crossing from one side of the surface to the other. 

The vortex strength (y> may be constant, as reported by Woodward and used in 
FLKXSTAh ( NAS2-5006), or have a more exotic variation such as the quadratic 
variation of the vortex spline (NAS2-6530). The vortex singularity has application to 
lifting problems e g , subsonic and supersonic and for steady ami nonsteady 

The vortex spline was introduced under NASA-Ames contract NAS2-H530. Pane Is uro 
defined over which quadratic strength variations exist in the spanwise direction and 
linear variations in the chordwise direction. A given spline function, defined by one free 
parameter governing its value spans a specified grid of neighboring panels such that the 
following conditions are satisfied for a given vortex spline function: 

• The value of the strength and slope (spanwise derivative) of the strength are zero 
on the edge of the grid except for possible special regions like the tips and the root 

• No discontinuities of value or slope occur w ithin the grid in the spanwise direction. 


• Strcumwise variations of strength an* linear it ml continuous 


Because <»f these conditions. although an inherent wake is present. there are no discrete 
trailing elements which produce infinite perturbation velocities in their vicinity In 
addition, perturbation velocities are continuous throughout 

A leust Minorca formulation of the boundary conditions wait used to obtain the free 
paramctcni characterizing the singularity strength variations The excellent success ol 
the method has previously been reported (ref 15). The method has been extended to the 
nonsteady subsonic flow regime under Task IV of NAS2-7729 with excellent results 

The use of any of the foregoing singularity types are conceptually identical In general, 
a singularity type is chosen for its ability to simulate a physical problem (see sec H O, 
Aerodynamic Modeling ». Having selected the singularity type, the values of the 
parameters governing its predetermined type of strength variation are established by 
imposing the boundary conditions on the boundary surfaces The process is conceptually 
independent of the type of singularity variation under consideration H e., constant 
strength, quadratic, etc ). 

The singularity methods being considered us candidate schemes include: (li FLKXSTAB 
aerodynamics (the reference point), (2) vortex spline (NAS2-6H30), 1 3 > source spline 
(NAS2-7729I, < 4 > doublet spline (NAH2-7729), (5) linearly varying vortex punels of 
Woodv/ard irefs 27 and 41 HH» Doublet Lattice (ref. 39), (7> line singularity source (refs 
I and 10). and (8) line singularity doublets irefs I and 11). 

Loading function methods such as that of Rowe trefs 1H and 42), although singularity 
methods, are of a type requiring the functions to span great regions of a geometrically 
restricted configuration and are not considered in this report These methods can, 
however, provide valuable comparisons for the methods considered here. 

7.2 KEATURKS OF CANDIDATE SCHEMES 

Some of the features of the candidate schemes will now be discussed Of interest is the 
formulation, the application to the boundary value problem, any restrictions imposed, 
computation features, and the future growth potential for the method The level to 
which the scheme hus b *en developed is discussed. 

Formulation 

The FLEXSTAB steady stute problem is form , ated using the source and vortex 
fundamental solutions, equations (H) and (Hi, with <r and y constant over ftanehi The 
source and vortex solutions share common paneling. 

The vortex panel has an associated wake which requires quadrilateral or triangular 
panels having side edges (possibly of zero length i in the direction of the free stream. 
This is illustrated in figure 7. 
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Figure 7 Constant Vortex Panel Geometry FL EXSTAB 


The panel geometry is defined in .« panel axis system such that the panel is in the plane 
Z 0. The houndury conditions have been linearised, all being expressed in the plane Z 
0. There is one boundary point per panel. 

The wake extends to downstream infinity oarallel to the X axis When a discontinuity 
of vortex strength or its span wise derivative is present, as in this formulation, infinite 
perturbation velocities occur in the vicinity of the discrete trailing elements (see fig. 7». 

Mathematically the vortex or pressure panel is a surface ncr • which a discontinuity ol 
tangential perturbation velocity exists The perturbation velocities are related to the 
vortex singularity strengths The perturbation velocities are applied to the boundary 
condition equations from which the unknown singularity strengths are determined 

The vortex spline scheme is built upon a similar foundation. The paneling is identical to 
thut of figure 7 However, instead of the singularity existing as a distribution over a 
single punel, each vortex spline, character lied by one free parameter, spans a grid of 
neighboring panels, generally a grid of 3 or 4 panels in spun and 2 panels in the 
stream wise direction. The spanwise variation is quadratic, the streamwise is linear (see 
figs H and 9l. 

A wake region exists in the vortex spline formulation but no discrete trailing elements 
are present because of the absence of jumps in singularity value or spanwise slope, i e.; 
no infinite perturbation velocities occur in the wake. 

The linearly varying vortex panels of Woodward (ref 27 > have a variation which is 
linear chordwise and constant spanwise. Kach of the singularity elements is 
characterized by on*- free parameter. Unlike the vortex spline, this method does have 
discrete trailing elements because of the discontinuous istepi spanwise variation and. 
therefore, infinite perturbation velocities do occur in the wake just as they do for the 
constant pressure panel method. A similar procedure is reported by Lopez und Sheri 
(ref. 43 1. 
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Figure 8 ■ Surface Spline Distribution of l orticity 



Figure 9 Spline Functions Distributed on a Wing I Symmetric Flow) 



The Houri »pllne ami doublet uplinv have no inherent wake Accordingly, the geometry 
requirement* imposed hv tin- vortex method* ure not necessary The only requirement i* 
that the panel* he quadrilateral or triangulur The source and doublet iqdinc* are 
therefore directly amenable to non linearised boundary condition* and arbitrary panel 
configuration*, in contra*! to the vortex method*. 

Mathematically the tuiurce and doublet panel* are. re*pectively, surface* aero** which 
either the perturbation normal velocity or the perturbation tangentiul velocity ha* a 
lump discontinuity Most often the source i* u*ed to *imulute thickne** effects and 
doublet* the lilting effects, however; for thick <nonlineari*ed> boundary surfaces, both 
type* together can u*uall> he used on the surface and or in the interior <>l bodies 

The ... /mutation of the source spline and doublet spline were presented in detail in the 
proposal submitted to perform this work 1 1 unit I Splined Venn>n of FI.KXSTAB 
Boeing document I)ii-ll7»>l April Itt74. submitted in response to Task I of HFP letter 
land February (• 1974. NASA Ames i. I’he source spline formulation ha* a linear 

variation in each of two directions Ixicul surface curvature is also included The Hourre 
strength* are not restricted to he continuous in value or slope at panel edges, but 
continuity can In- very nearly enforced hv the solution. 

For subsonic llow, the three-dimen*ional doublet spline formulation* lor analysis and 
design ar< not strongly geometry dependent. The doublet distribution* are continuous 
quadratic functions over panel* hut (from panel to panel i are not necessarily continuous 
in vuluc and slope at ail points of the panel edges The doublet spline formulation leads 
to linear equation* tn the unknown doublet strengths which enables n linear influence 
coefficient formulation of the entire problem Surface curvature effect* are included 
The doublet spline formulation has led to a simplified form for the influence 
coefficient* Because the doublet distribution ir represented by a simple polynomial 
(quadratic surface function i. all integral expression* for the influence coefficient* can be 
integrated in dosed form and with considerably fewer terms. Thi* produces less 
computation time and provide* greater numerical reliability. 

For supersonu flow, supersonic doublet Splines are being developed as > part of the 
Task III ( NASH-77291 studn - Studies of the supersonic doublet splines were earlier 
initiated under NASA-Ame* contract NAS. H530 In that study, the doublet 
characteristic box method was introduced Briefly, the method has the following 
feature* First, special Mach lines emanating from planforin edge break* are identified 
These lines ore used to divide the plunform into several different region* and each of 
these region* i* divided into a network of panel* hounded by Mach lines. Kach panel 
ma> then fall into a relatively small i 10(, restricted class of basic numerical building 
block element*. A planforin i* represented by a superposition of these elements. A 
quadratic doublet singularity distribution, having a specific number of free parameters, 
is assumed to exist over each panel or building block element. The unknown free 
parameters are determined hv i .lunations of downwash conditions and vorticity 
continuity conditions Application of tin* method to a cambered triangular region 
hounded by a supersonic leading edge and one *peciul Much line ha* produced results 
which are indistinguishable from exact theory. 


Thi* doublet lattice (fb; 1 0* im a method applicable only to subsonic (low A line of 
acceleration potent ml doublet* lie* along a panel quarter chord and ha* an unknown 
< constant > singularity strength One rolloeafion point, located at the panel mid-*paii ami 
at the three-quarter chord, m used to determine the strength; this boundary point 
position generally differ* from •*»««» unit! for the mmlifting problem Trailing vortex 
element* are included to allow for lifting solutions The doublet lattice method ha* 
direct application to the arbitrary frequency nonsteady problem iref 39» and ha* been 
applied to a wide rarge ol appli minus le g., ref 44 > 



The line source und line doublet method* distribute source* and doublet* along an uxial 
line, the variation of which may be constant, linear, quadratic, etc , in a coordinate 
system representative nf the line Boundary points are located at selected points along 
the line at t!.e real (or imagined) surface. Figure 1 1 illustrate* the line singularity 


formulation. 
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Figure 1 1. ■ /.»/»* Singularity Formulation 

Boundary Value Problem Application 


The vortex methods ar< generally used with linearized boundary conditions because of 
the inherent wake The source and doublet methods are used with linearized boundary 
conditions and also with exact boundary conditions in subsonic (low. For supersonic 
(low. Woodward (ref. 27) has been successful in using sources on a body surface with 
exact boundary conditions. The work of Task III (NAS2-7729) has investigated source 
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and doublet spline surface paneling in supersonic flow The lint- singularity method* 
employ I neariwd boundary conditiomi. 

The vorte- methods ure more restricted in paneling requirements than the source and 
doublet Of the two candidate vortex methods, FLKXSTAB's constant pressure panels 
have the most severe paneling restrictions, having not only the wake imposed 
restrictions on downstreum paneling but also restrictions on the relative sire of 
neighboring panels and the location of the boundary points. The vortex spline does not 
suffer as markedly from any of these paneling restrictions. However, it may be more 
difficult to set up the geometry of a general configuration for the vortex spline. The 

restrictions seem to be of minor concern for the source and doublet me hods Cases have < 

been run using the doublet spline for which panel corner points were generated in a 
random manner with no loss in computational accuracy compared to cases run with 
more conventional paneling 

The use of doublets to simulate a luting configuration would require a paneled wake 
geometry across which a jump in potential (doublet strength) can occur. This is an 
advantage when the geometry of the wake is important as in high lift configuration 
applications Only u few panels may be necessary to represent the entire wake, e g , in 
two-dimensional flow, one semi nite doublet panel can represent the entire inhe.ent 
planar wake produced by the vorte,* method 

A note of caution is in order concerning the vortex spline method If the geometric panel 
grid is planar, the vortex spline produces no discrete trailing elements. If the geometry 
has curvature, the curvature must be continuous in value and span wise slope to avoid 
discrete trailing elements produced by the geometric discontinuity This not only 
imposes additional geometry paneling requirements, but also generates additional 
complexity within the influence coefficient kernel functions (see eq. (8)1 Whereas the 
other candidate schemes require only single boundary points per panel, the vortex 
spline may use the method of least squares which requires jur boundary points per 
panel in supersonic flow and two in subsonic flow. This increases the size of the 
influence coefficient matrices and, accordingly, the computation of either a matrix 
solution or a matrix inversion, ore of which is necessary for the analysis problem. Since 
the constant pressure panel method of FLEXSTAB has the discrete trailing elements in 
general, whether they are produced by discontinuities in singularity strength or 
geometry, and it, therefore, can be applied to configurations of arbitrary spanwise 
geometry when proper account is taken of the paneling requirements 

Any of the singularity schemes presented may be applied to the solution of the 
appropriate nonsteady boundary value problem The governing diffeiential equation and 
boundary condition are given by the linear equations (ll and (2). A linear relationship 
can be found relating the zero-order potential solution (steady-state flow) and the 
first-order in time potential through a nonhomofeneous linear boundary value problem. 

Computation Features 

A quantitative comparison of the computation features of the candidate methods is 
difficult. Certainly the efficiency of the coded algorithms determines the length of the 
computation time for any of the methods In addition, the form of the integration (i.e.. 
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unHly'.ial and nr numcrnal quadrature of the integral expressions), even if efficiently 
coded, mas not prove to be the fastest manner of computation possible Table '2 has been 
assembled for the purpose of illuslratir, the items which alfect the computational speed 
and efficiency alonK with other features of the various schemes Note that some items, 
such as surface curvature, may increase a local computation cost but may decrease (he 
total * global • cost because less f jnels are required 


Table 2 Comparison Features of Candidate Numerical Schemes 
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Current Level of Development (August IH74) 

The FLKXSTAH constant vortex panel aerodynamic influence coefficient p-ogram in 
fully operational Thin includcH routines to perform the low frequency nonstnady 
analysis for noriplanar configurations The vortex spline has been formulated for the 
steady flow, planar configuration case with some work in subsonic nonsteady flow. In 
subsonic flow, both the source and doublet splines are formulated for steady nonplauar 
flow In steady supersonic flow, the doublet and source spline are under development 
The doublet lutlice method isubsonic flow only) is well developed and documented for 
steady and imnstcady flow The line singularity method is similarly available for both 
subsonic and supersonic flow. 

Future Growth Potential 

The constant pressure panel has been successfully used for the subsonic and supersonic 
flow regimes, both steady and low frequency nonsteady. For steady flow, a linear 
chordwise variation has been introduced <ref. 27 1. The vortex spline is considered less 
general thun the constant pressure panel due to the configuration restrictions 
previously discussed. 

The source and doublet spline methods ure proving to be extremely powerful and have 
application to a wide range of problems To illustrate the level at which the Bource and 
doublet splines are performing, it was the numerical building block used in the work for 
NASA -Langley contract NAS1-12185 and its follow-on contract NAS1-13H33. The 
objective was to conduct an analytical study, develop a method of solution, and develop 
a computer program to predict the subsonic aerodynamic loads on a delta-like wing over 
which a leading-edge vortex exists. This is a very complex nonli»ear mixed analysis and 
design problem requiring an iteration procedure for solution. The source and doublet 
spline successfully hundle this problem, whereas the vortex panel methods have little 
chance for success. The source spline with curvature has recently been used to calculate 
the subsonic potential flow over a sphere. A constant source panel method (ref. 32) 
required 1342 elements while the source spline required only 1H2 elements to achieve 
the same accuracy. 

Figure 12 is a qualitative illustration of where the various panel singularity methods 
fall on a line of computational accuracy. Note that some items on the lius are not 
attainable as, for example, exact solutions are difficult to obtain (one of the only known 
exact three-dimensional solutions is for incompressible flow over a circular wing). 
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Figure 12. -Panel Singularity Methods 
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H.o AKKODYNAMK MODULI Nfi/NUMKRK'AL SCHEME 

COMBINATIONS 


It in the purpose of thin section to make recommendationH from among the candidate 
schemes of sections H.O and 7.0 and to assess the impact of the recommended 
combinations on the ultimate objectives of Task I An estimate is made of the resources 
required to achieve those objectives. 

H I RECOMMENDATIONS 

Before recommendations can he made, it is necessary to consider the particular version 
of FLEXSTAB to which they relate. This was briefly introduced in section 4 2 The 
present FLEXSTAB system is an intimate collection of programs. A revision of any 
magnitude generally impacts the entire program. This is because the original program 
was formulated without any particular regard to a modularized structure The condition 
was precipitated bv taking off-the-shelf programs (primarily the Alt' programs! and 
forcing conformity at another level of the FLEXSTAB system, for example, the 
aerodynamic programs and internal structural influence coefficient and externul 
structural influence coefficient programs (IS1C and ESICi were interfaced within the 
stability derivatives and static stability programs. This has produced a close-coupled 
system of programs for which the general interface problem was never fully addressed. 
A change within the aerodynamics program is sure to affect portions of the SD&SS 
program This formulation of FLEXSTAB is here labeled a Level 1 version for later 
discussion. 

A • *>vel 2 version of FLEXSTAB is one which has h..d careful attention puid to a 
modular structure type of formulation In simplistic form, one routine could be 
unplugged and another substituted in its place. A Level 2 version of FLEXSTAB would 
not have the inheunt weaknesses of the Level 1 version. It would be able to 
dynamically participate in the changing environment of advancing technology w ith the 
flexibility to specialize it for specific applications. TF interfaces at which fundamental 
duta flow takes place is of primary importance for Le.el 2 work This data flow should 
be clearly identified and the particular format specified. Level 2 provides the best 
chance that advancing technology can be received into the FLEXSTAB en’ i<onment 

A generul flight vehicle configuration is con. posed of many items: wing surface, body, 
control surface, nacelle (engine), strut, wake and or exhaust region. It is considered 
advantageous to have various degrees of component modeling. This enables flexibility of 
application for the users and a gradation of accuracy for the various levels of analysis 
common to the aeroelastic design cycle, from preliminary design through final 
configuration analysis. 

For purposes of discussion, a notation has been developed to allow easy reference to the 
various candidate aerodynamic (table 3) and numerical schemes (table 4) presented in 
the preceding two sections. Only certain combinations of these aerodynamic models and 
numerical schemes are possible. Table 5 illustrates them with an X. Within table 5 is 
shown the present FLEXSTAB combination for wing-body analysis <B.1:LS. Lll; W 1 


SS (constant), VC; 1 1; VC; N I VC). The body m represented by line Hingularity 
distributions of line sources and line doublets. The wing is a mean surface wing of 
constant source and vortex panel, and the interference shell is composed of constant 
vortex panels, as is the flow through nacelle 


Table 3 Aerodynamic Modeling Schemes 


• Body 

B. 1 

Lme representation tor isolated body (♦ 1.1, and 1.2) 

B.2 

Surface panel representation (+ 1.3) 

• Wing 

W.1 

Mean surface 

W.2 

Exact surface 

• Intel feience 

1.1 

Intel terence shell of constant CS 

1.2 

Interference shell of nonconstant CS 

1.3 

In'ernal lifting surface 

• Nacelle 

N.1 

Flow through (either W 1, W 2, 1.1, or 12) 

N.2 

Streamline nacelle (appropriate B.1) 

N.3 

Engine on (W.1 W.2, 1.1, oi 1.2) 

• Exhaust and wake 

E.1 

Specified wake 

E.2 

Jet exhaust wake 


Table 4 Numerical Schemes 


e Line 


LS 

Line source 

LD 

Line doublet 

• Source Panel 


SS 

Source spline 

• doublet panel 


DS 

Doublet spline 

DL 

Doublet Lattice 

e Vortex panel 


VC 

Constant 

VS 

Vortex spline 

VL 

Linear m-chord 
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Table b Possible Combinations of Numerical Schemes and Aerodynamic Models 



Line 

Source 

Doublet 

— 

Vorte* 


LS 

LD 

SS 

OS 

Dl 

VC 

Vo 

VL 


B. 1 

X 

X 







Body 











B.2 



X 

X 

X 





W.1 



X 

X 

■TH 


X 

HE'S 

Wing 

W 2 



X 

X 

M 

hi 


m 


U 




X 

X 

X 

X 

X 

In ter fur- 

1.2 




X 

X 




ence 

1.3 




* 

* 


Planar only 



N.l 



X 


Linearised boundary condition only 

Nacelle 


X 

X 


s 









X 


X 

X 

X 

X 

Wake/ 

El 





X 


Planar only 


exhaust 

E.2 



X 



Constant cross section 



W 1 



X 

mm 

mm 

X 

X 

X 


B.2 



X 







♦ 










1.3 







Planar only 













W 1 



X 

X 

X 

X 

X 

X 


♦ 










B. 1 

X 

X 








1.1 



X 

X 

X 

X 

X 

X 


Of 










1.2 




X 

X 

X 

X 

X 


Two option* art- now presented which are aimed at the two particular version level* 
previously discussed lor the KI.KXSTAB system. Within each option, recommendation* 
are presented to achieve the ultimate objectives of Task 1 Within the recommendation* 
there is riot yet a clear choice of one numerical scheme over anc-iier. There are four 
reasons for this 

1. Net all candidate schemes (e g., supersonic doublet spline) are sufficiently 
developed to warrant an unqualified acceptance or rejection. Such schemes are 
being included with the idea that they will be more fully evaluated at a later time. 

2. Different levels of modeling combination* and levels of accuracy should be 
available to a user to more adequately match his immediate need and the 
knowledge he has of his configuration. This mav be possible under option B. 

3. Verification of the recommended combinations is discussed in section 8.3. The 
studies conducted at the time of verification will not necessarily eliminate a 
particular modeling combination, but instead will serve to establish the impor'ant 
comparison characteristics of the several modeling combinations. 
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4 ' r ht> impact of the modeling combination* on the low frequency capability is 
discussed in section H.2 in a qualitative manner. Such discussion can only be 
quantitatively established at the time of verification (see point 3). 

Because the formulation* of the doublet uplines and source splines include higher order 
singularity strength variations and surface curvature and because the other schemes 
are subsets of these, they are to be considered the preferred schemes. 

Option A 

Option A assumes that only a Level 1 version of the Fl.KXSTAB system is available, 
that is, the present Fl.KXSTAB will not undergo a restructuring procedure 
Accordingly, the entire aerodynamic modeling of Fl.KXSTAB is retained, i.e.; line body, 
mean surface wing, interference shell, etc The use of an exact body aerodynamic model 
in the Level I version is considered unwise. The numerical methods recommended under 
option A are the new technology methods emerging from the work of Task 111 
(NAS2-7729) along with the vortex spline ( NAS2-6530). These numerical methods would 
he applied to the wings and shells. 

Several key points can be noted regarding the numerical methods. Source splines are 
used to simulate thickness effects for the (linearized boundary condition) wing. Vortex 
splines or doublet splines may be used for lifting effects of the wings and shells, 
although the choice of doublet splines allows the flexibility to use the specified 
nonplanar wake model (K 1) The li near-in-chord vortex method is not recommended 
because it is not a significant improvement over the constant pressure panel method, 
and is a subset of the vortex spline The recommended combinations for option A are 
shown in lublett 


Tdhh> 6. Option A Recommendation of Combinations for Wing Body Analysis 



Line 

Source 

Doublet 

Vortex 

LS 

LD 

SS 

DS 

DL 

VC 

VS 

B.1 

X 

X 






W.1 



X 

X 

X 

X 

X 

1.1 




X 


X 

X 

1.2 




X 





Option B 

Option B assumes that Fl.KXSTAB will be restructured to the status of a Level 2 
system of programs The exact body surface aerodynamic model is recommended <B 2) to 
provide an improved modeling capability not presently available. The selection of the 
exact body model all but eliminates the vortex numerical methods for the exact body 
representation because of their inherent wake, but does not necessarily eliminate them 
for other modeling applications (e.g., wing representation). 
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The aerodynamic model recommended for the wing is the mean surface representation 
presented with option A and lined within the present FLEXSTAH system This modrl 
employ* linearized boundary conditions and include** thickness and lifting effects. Thin 
analysis ia generally adequate f.*r sections up to 15** thick For aeroelastic analysis, the 
load ih important. Surface pressures may also he obtained 

Apart from the direct application for aeroelastic analysis, there is another point which 
should be considered This •» the calculation of lateral-directional derivatives (ref. 45i, 
When the wing dihedral is very near zero und the linear representation is used, it is 
dill u ult to compute certain 'oupled stability derivatives for nonzero wing dihedral, the 
method shows no su<'h difficulty and those particular stability derivatives are ol 
second-order 

The exact surface representation ol the wing (W.2» would exhibit no difficulty for near 
zero dihedral unlike the linearized one. However, it has not been firmly established 
that the linearized model could not be properly formulated to remove this difficulty. 
This should be investigated. 

In addition to its direct application to aeroelastic analysis, the mean surface 
representation would result in smaller matrix sizes and therefore faster computation 
times 

The numerical methods recommended for the mean surface wing are the source splines 
for thickness and vortex nr doublet splines for the lifting effects. The use of doublet 
splines allows the flexibility to use the specified nonplunur wake model (El) and is 
thus preferred The linear-in-chord vortex method is not recommended because it is not 
u significant improvement over tin* constant pressure panel method, and is a subset of 
the vortex spline. The combinations for option 11 are show n in table 7. 


Table 7. Option B Recommendation of Combinations for V\/iny Body Analysis 



Line 

Source 

Doublet 

Vortex 


LS 

LD 

5S 

DS 

DL 

VC 

VS 

B 2 



X 





W.1 



X 





1.3 








Nacelle, Wake, and Exhaust Combinations 

Both options A and B can employ the following combinations within their analysis. A 
linearized boundary condition formulation is suggested for option A; either exact or 
linearized boundary conditions are suggested for option B. 

Three alternate nacelle aerodynamic models were presented in section H t): flow through, 
stream tube (entrainment), and engine-on For models N.l and N.3. the source spline is 
recommended for thickness effects and the vortex spline (linearized boundary 
conditions) for lifting effects. The use of the doublet spline ullows the flexibility to use 
the specified, nonplanar wake model (E.li. The nacelle model <N.2i employs a line 













source Httd thereby in able to simulate entrainment effor'i (a* an entrainment model it 
is more properly termed a jet exhaust model). These alternate m (dels allow a range of 
capability and modeling sophist nation not presently available in FLEXSTAB Nacelles 
no longer need to be circular nor to In- of constant cross section. 

The wuke aerodynamic model «E I* is specified by its geometric description, the planar 
wake being a subset of the nonplanar Because of the absence of an inherent planar 
truilmg wake the doublet methods are recommended for the numerical scheme 

The jet exhaust model iK 2> employs either the doublet splines and or source splines for 
its numerical description The plume shape can be specified by the user or set by a 
predetermined default option. 

The combinations just described ure shown schematically in table 8. 

Table 8 Recommendation of Combination for Nacelle Wake and Jet Exhaust Analysis 



Line 

Source 

Doublet 

Vortex 


LS 

LD 

SS 

DS 

DL 

VC 

VS 

N 1 



X 

X 

X 

X 

X 

N.2 

X 







N.3 



X 

X 

X 

X 

X 

E.1 




X 




E.2 



X 

X 





H.2 DISCUSSION OF TASK I OBJECTIVES 

Each of the combinations recommended in section 8.1 (especially option Hi. adds 
markedly to the present aerodynamics capability of the FLEXSTAB system. Because the 
constant pressure panel has not been eliminated from consideration it is clear that no 
capability, relative to the NASA Ames released version 1.01.00, has been lost The 
combinations offer additional variety for modeling and accuracy levels not yet available 
in the FI, EXSTAB system 

The present nonstendy capability w.thin FLEXSTAB is limited to the low-frequency 
approximation None of the above combinations is expected to cause a loss ot th«s 
capability. In fact, the studies of Task IV < NAS2-7729) with the vortex spline have 
demonstrated conclusively its applicability to the full nonsteady problem The doublet 
methods will no doubt also have application to the more general arbitrary frequency 
and Kussner- Wagner formulations of the nonsteady problem (ref. 46). This of course will 
need to be verified outside of the present study for the FLEXSTAB environment 

It is possible to formulate the aerodynamics program in a manner to guarantee that it 
can be interfaced with other program types (e g., structures) This will require much 
work to insure that appropriate interface data is made available for the wide range it 
potential uses The present FLEXSTAB is not formulated with the functional module 
concept, the aerodynamics, structures, geometry, etc., are intimately related This point 
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Mill do ter mint* the future course ol tlx- present study and the ultimate fulfillment of the 
Tank I objectives It m outside the scope of the present study to offer recommendations 
as to how FLKXSTAB should he structured Unless a major reorganization is performed 
to formulate the functional module con 'opt lor FI.KXSTAB (Level 2 version i, little 
chance is given that FI.KXSTAB Mill m'cruct dynamically with advancing technology 

Retention ol the Level I version of the program will dictate that any improved 
aerodynamic capability will have to la- converted into an equivalent constant pressure 
panel formulation at the local aerodynamic structures interface The only payoff from 
an improved aerodynamic program in a Level I program will he the output data from 
the aerodynamic program, with little improvement to be observed downstream ol that 
program It is virtually impossible to incorporate a second order lateral directional 
capability into the Level I program. 

The option A payoff is a minimum impact on the FI.KXSTAB system and low technical 
risk Bern u •* o| the minimum impact, a minimum time would be required to develop a 
working program assuming the existence of the low frequency formulation and 
presently nonexistent routines such as the supersonic doublet numerical scheme 
Numerically, there would Ik- an expected increase in accuracy, elimination of control 
point sensitivity, elimination ol paneling sensitivity, and a reduced sensitivity to tail 
arrungenn ni In addition, a working version of the FLKXSTAB system, which 
incorporates the improved numerical techniques, would In- available quickly und with a 
minimum ol effort. The limite.tions of option A include the retention of the crude 
aerodynamic model ( interference shell, etc ! which brings into question the advisability 
incorporating improved numerical methods and the relative impossibility of a 
:mic program capable of accepting advancing aerodynamic technology. 

The significant payoff from option B is the ability to model the actuul configuration, a 
point considered extremely important for quality results from a potential flow program. 
Associated w ith this development is a restructuring of the FLKXSTAB syste n the 
effect of which is to impact heavily the time required to obtain a new ver> ion of 
FLKXSTAB The corresponding technical risk is higher, not in terms of technical 
feasibility but in terms of development time because unforeseen difficulties often cause 
schedule slides. Accordingly, the investigation und development of the restructuring 
procedure und development ol individual modules must be done with due regard to 
careful, parallel development and with special attention to seemingly infinitesimal 
detail The anticipated numerical improvements hold the promise of an aeroelustic 
capability of a quality never before achieved. 

H..’f VERIFICATION PROCEDURES 

Verification is intended to provide quantitative grounds for comparison of the 
combinations and to define their respective regions of application Verification will not 
generally establish superiority of one combination over another in all characteristics 
(e g., accuracy, efficiency of computation, absence of geometric constraints, applicability 
to general configurations, etc.i. Instead, verification supplies confirmation that the 
modeling combination (aerodynamic model and numerical schemei can supply the key 
features required to model the physical flow problem of interest to the user Verification 
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ih not a part of thin contract hut would hi* one of the item* of work to help fulfill the 
funk I ultimate objectives 

Slumlord* of Comparison 

Three standards of comparison for a combination are: exact analytic theory, other 
existing numerical methods, and experimental data The nature of experimental data 
often makes it a good qualitative comparison lor a potential flow solution. Kxact 
solutions of the boundary vulue problems would he the ideal standard of comparison hut 
these are not generally avniluhle. 

Of the three standards, comparison of the combinations with other existing numerical 
methods oiler the most flexible and valid verification for those combinations. Such 
methods do exist and their validity is widely recog nir d These include the methods 
reported in references Ifi, 16, 23, 26, and 27 

Variables of Comparison 

Many quantities can he compared. The pressure distribution on a configuration, 
however, is the primary aerodynamic quantity of interest, Other quantities such as 
total lift and moment, sectional lift and moment, etc., are integrated quantities ana 
may mask difficulties inherent in a particular scheme. These should not he considered 
primary quantities Data should also he compiled demonstrating the economics of using 
a particular combination. This includes computation time to achieve a specified level of 
accuracy. 

Items to be Verified 

The following items should be addressed as part of the verification procedures for the 
various combinations 

• What are the basic spline formulations required in regions such as wing-body 
intersection, w ing planform breaks, and wing tips? 

• Do the combinations have application to the low frequency nonsteady flow 

problem? 

• Does the combination exhibit numerical stability? 

• What is the sensitivity to panel size, arrangement, and control point placement? 

• Does the solution converge for increasing panel density? 

• Are the combinations adequate to model the physical flow? 
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hollow -On Item* of Work 


Hu follow -on itcniM ol work associated with option* A and If ure 

Option A 

1 Make « comprehensive evaluation un<i verification against the preaent KLKXSTAB 
and against more refined at .(dynamic model* <e g . ref. 23). 

2 Add re** tin* KLKXSTAB interface problem and the converninn to a constant 
pressure panel equivalent. 

3. Make decision to proceed or scrap 

Option It 

I Validate the aerodynamic model* and numerical method* recommended. 

2. Address the KLKXSTAB interface (must be done in parallel with any KLKXSTAB 
restructuring plan*) 

3. Make decision to proceed or scrap 

8.4 KKKOI KCKS KKQUIKKD 

The estimated resources (August 1P"4 figure*) required to achieve the ‘ndividual items 
ol work toward the fulfillment of the ultimate objective* of Task I are presented below 
I'he remaining work is broken into phases (2 through 4) and expressed individually for 
each option )A and B) Note t hint for option B. no cMtimute is made of the resource* 
required to develop a Level 2 version of KLKXSTAB All resource estimates are 
contained in table t) 

Option A 

The assumption* made to generate the figures contained in table t» for option A are 

• Source splines and doublet splines will be the fundamental numerical building 
blocks. 

# The subroutines will he available for steadv stuti doublet und source splines, 
subsonic and supersonic 

* The low frequency doublet building blocks will bo available 

s A converter will be developed to obtain constant pressure panel equivalents. 

* The present KLKXSTAB aerodynamic modeling is retained 
The manpower level is 4 . r iH manyeurs 
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Table 9 Resour c et Required to Fulfill the Ultimate Objectives of Task / 


August 1974 figuret 





Manmontht 

Manmontht 

Computer ret. 

Flow tim« 




BCAC a 

BCS b 

CU’l 

montht 




Verify 

8 


4 


200 


8 




2 

Interface 

2.6 


26 




6 

8 




Modeling 

4 


06 


- 


4 



Option A 

3 

Converter 

3 

30.6 

0.6 

24 5 

- 

434 

3 

4 

25 


Program 

3 


10 


20 


10 




4 

implement 

2 


2 


4 


2 

13 




Verity 

4 


1 


200 


4 




Document 

6 


4 


10 


5 





Verify 

10 


6 


200 


10 




2 

Interface 

2 5 


25 




5 

10 




Modeling 

6 


0 5 




6 



Option B 

3 

Interface 

2 

34.5 

1 

30 

— 

454 

2 

6 

34 



Program 

3 


12 


40 


12 





Implement 

1 


1 


4 


1 

18 



4 

Verity 

4 


1 


200 


4 





Document 



6 


10 


6 




a BCAC Boeing Commercial Airplane Company 
1> BCS Boeing Computer Services, Inc. 


Option H 

The- assumptions made to generate the figures contained in table 9 for option B are: 

• Option B costs are independent of a parallel, Level 2 development. 

• Source and doublet splines will be the fundamental numerical building blocks. 

• The su w - lutines will be available for steady-state doublet and source splines, 
subsonic and supersonic. 

• The low frequency doublet building blocks will be available. 

• The exact surface representation will be the aerodynamic model for the body. 

The manpower level is 5.38 manyears. 


PAGE W 


9.0 CONCLUSIONS 


The ultimate objectives of Tank I can be achieved with little technical rink and without 
Iok* of any capability relative to the NASA Ame* released version 1.01.00. The direction 
of that development and the ability of the FLEXSTAB computer program system to 
interact dynamically with aovancing technology hinges upon development of a 
restructured program system. Option B is recommended if the restructuring is done; 
otherwise, option A is recommended, the payoff versus the cost of which is questionable. 
Option B is the specific recommendation of this report. 


Boeing Commercial Airplane Company 
P.0 Box :»707 

Seattle, Washington 9H124 
August 31, 1974 
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